NASA Technical Memorandum 106535 
ICOMP-94-04 


/ 


Theoretical Study of the Incompressible 
Navier-Stokes Equations by the 
Least-Squares Method s 

O 

O 

I 

O' 

z 


rg 

m O 

«T> O 

r— m 

U O 

C O 

3 O 


m 

m 

O 


Bo-nan Jiang 

Institute for Computational Mechanics in Propulsion 
Lewis Research Center 
Cleveland, Ohio 


Ching Y. Loh and Louis A. Povinelli 
National Aeronautics and Space Administration 


> 00 
Q Lu 

3 

o 

O0 ►— 

00 00 

J I lAJ £ 

< cc a: u 

UUJ<U 
M H 3 ft! 

k > a a> 

uu < 00 in 
at 2 | fU 
O >-* a: 

LU UJ OO 

X -J < tn 
K CQ Hi 

-J 3C 
00 a; 

^ LO Uj U 
lAUJ I 


Lewis’ Research Center 
Cleveland, Ohio 


March 1994 


man- • 

in CL < CL 
%n >: v 00 
O O CO <£ so 

rH O Z s* 

| ? LO ^ 

X z 

►- a ^ 

| UJ »"• o u 

< x o a) 

to h < I 

< 3 C 

ZHOUi li 
w O UJ X o 



National Aeronautics and 
Space Administration 




Theoretical Study of the Incompressible Navier-Stokes Equations by the 

Least-Squares Method 


Bo-nan Jiang 

Institute for Computational Mechanics in Propulsion 
Lewis Research Center 
Cleveland, Ohio 44135 


Ching Y. Loh* and Louis A. Povinelli 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


Abstract 

Usually the theoretical analysis of the Navier-Stokes equations is conducted via the 
Galerkin method which leads to difficult saddle-point problems. This paper demonstrates 
that the least-squares method is a useful alternative tool for the theoretical study of partial 
differential equations since it leads to minimization problems which can often be treate 
by an elementary technique. The principal part of the Navier-Stokes equations in the 
first-order velocity-pressure- vorticity formulation consists of two div-curl systems, so the 
three-dimensional div-curl system is thoroughly studied at first. By introducing a dummy 
variable and by using the least-squares method, this paper shows that the div-curl system 
is properly determined and elliptic, and has a unique solution. The same technique then is 
employed to prove that the Stokes equations are properly determined and elliptic, and that 
four boundary conditions on a fixed boundary are required for three-dimensional problems. 
This paper also shows that under four combinations of non-standard boundary conditions 
the solution of the Stokes equations is unique. This paper emphasizes the application 
of the least-squares method and the div-curl method to derive a high-order version o 
differential equations and additional boundary conditions. In this paper an elementary 
method (integration by parts) is used to prove Friedrichs’ inequalities related to the div 
and curl operators which play an essential role in the analysis. 


♦National Research Council-NASA Research Associate at Lewis Research Center. 



Introduction 


The incompressible Navier-Stokes(NS) equations which govern the motion of viscous 
fluids are among the most important partial differential equations in mathematical physics. 
The Navier- Stokes equations are derived from the mass and momentum conservation laws 
and the linear local stress-strain relation. While the physical model leading to the Navier- 
Stokes equations is simple, the situation is quite different from the mathematical point 
of view, and can be described by quoting the words from a mathematician: “Despite the 
important work done on these equations, our understanding of them remains fundamen- 
tally incomplete.”(see Teman[l]). As pointed out by Teman, the major difficulty in the 
mathematical study of these equations is related to their nonlinearity. 

Even our understanding of the linearized Navier-Stokes equations and the Stokes equa- 
tions remains incomplete. Here we refer to the fact that there has been no systematic study 
of permissible boundary conditions(BCs). The boundary conditions applied to the Navier- 
Stokes equations have been the subject of constant controversy. Quoting next from a 
recent review on incompressible flow, Gresho[2] states that “To the best of our knowledge, 
the jury is still out regarding the full story on mathematically permissible BCs for the NS 
equations” . 

Looking into the classic books on the mathematical theory of the Navier-Stokes equa- 
tions, see for example, Ladyzhenskaya[3] , Teman [1,4], and Girault and Raviart[5], one will 
find that mathematicians study the existence and uniqueness of the solutions mainly for 
the standard case in which the velocity components are prescribed on the boundary. On 
the other side, physicists and engineers apply the non-standard boundary conditions often 
based on physical intuition. 

Usually the Galerkin method is employed for numerical solution and theoretical analy- 
sis of the incompressible Navier-Stokes equations. The application of the Galerkin method 
leads to a saddle-point problem, and thus a difficult “Ladyzhenskaya-Babuska-Brezzi(LBB) 
condition” or “inf-sup condition” is involved [6-10]. 

In recent years the least-squares finite element method(LSFEM) has been success- 
fully developed and used for the computation of incompressible viscous flows[ll-19]. The 
LSFEM based on the first-order velocity-pressure-vorticity formulation has superior ad- 
vantages over other methods: together with, for example, a Newton linearization, the 
resulting discrete equations are symmetric and positive definite, and can be efficiently 
solved by matrix-free iterative methods; the choice of approximating spaces is not subject 
to the LBB condition, and thus equal-order interpolation with respect to a single grid for 
all dependent variables and test functions may be used; only simple algebraic boundary 
conditions are imposed, no artificial numerical boundary conditions for the vorticity need 
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be devised at boundaries at which the velocity is specified; accurate vorticity approxima- 
tions are obtained; neither upwinding nor adjustable parameters are needed; all variables 
are solved in a fully-coupled manner, operator splitting turns out to be unnecessary. In 
addition, we found that the least-squares method based on the first-order velocity-pressure- 
vorticity formulation provides not only a powerful technique for numerical solution but also 
a useful tool for theoretical analysis of the Navier-Stokes equations. 

In our previous work[l5], we showed that the original three-dimensional velocity- 
pressure- vorticity formulation for incompressible Navier-Stokes problems is not elliptic 
in the ordinary sense, and the compatibility condition, that is, the solenoidality of the 
vorticity, should be added to make the first-order system elliptic. Based on the first-order 
elliptic system we showed that for three-dimensional incompressible viscous flows four 
(three in most cases if the condition for the dummy variable is not counted) boundary 
conditions are required for a fixed boundary. In this paper, we give a further detailed 
discussion of the permissible non-standard boundary conditions, and prove the existence 
and uniqueness of the solution under these boundary conditions by using the least-squares 
method. The least-squares method leads to a minimization problem which is much easier 
than a saddle-point problem and often can be dealt with by an elementary technique. 

According to the explanation in [15], the non-linear convective terms in the Navier- 
Stokes equations have no effect on the classification. The ellipticity (in the ordinary sense) 
of the Navier-Stokes equations is determined only by the principal part of the equations. 
The principal part of the Navier-Stokes equations is the same as that of the linear Stokes 
equations. Also according to the theory in [20], error analysis for the nonlinear Navier- 
Stokes equations is essentially the same as that for the linear Stokes problem, at least 
away from singular points. Thus, the main goal of this paper is the verification of the 
well-posedness of the boundary conditions for the Stokes equations. 

The Stokes equations written in the first-order velocity-pressure-vorticity formulation 
consist of two coupled div-curl systems. The three-dimensional div-curl system is tradi- 
tionally considered as an overdetermined system, and its numerical solution is not trivial, 
since there are three unknowns and four equations. It is worthwhile to point out a similar 
situation in electromagnetics. Maxwell's equations also consist of two coupled div-curl 
systems. For the same reason, some engineering books on electromagnetics claim that 
Maxwell’s equations are overdetermined, and only two curl equations are independent; in 
some works on computational electromagnetics the divergence-free equation of the electric 
or magnetic field is just ignored. This paper demonstrates that it is incorrect by simply 
counting the number of unknowns and equations to judge whether a system of differential 
equations is overdetermined or not. Because of the importance of this problem, in the first 
part of this paper we show that the div-curl system is properly determined and elliptic, 
and analyse the div-curl system by using the least-squares method. In the second part of 
this paper we use the same technique to deal with the Stokes problem. 
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In fluid dynamics as well as in electromagnetics there are some good reasons why 
other higher-order versions of the Navier-Stokes equations and Maxwell’s equations are 
often useful (see the discussion by Gresho[2]). These are all derived from the “primitive” 
equations by differentiation and often offer additional insights regarding fluid flow and 
electromagnetic fields. They also serve as the starting point for devising alternative nu- 
merical schemes. For example, the pressure Poisson equation is obtained by applying the 
divergence operator to the momentum equation in the velocity-pressure formulation; the 
vorticity transport equation is obtained by applying the curl operator to the momentum 
equation. A key issue is that a derived equation obtained by simple differentiation admits 
more solutions than do its progenitors, in other words, spurious solutions (or spurious 
modes in electromagnetic waveguides) may be generated. With careful selection of addi- 
tional boundary conditions (and additional equations), the solution of the derived equations 
will also solve the original equations. In this paper, as a by-product, we show that the least- 
squares method provides a systematic and consistent way to derive higher-order versions of 
differential equations and corresponding boundary conditions without generating spurious 
solutions. That is, the Euler-Lagrange equations associated with the least-squares weak 
formulations are the most appropriate derived equations. We also show an alternative sys- 
tematic method to derive equivalent higher-order systems. This div-curl method is based 
on the theorem: if a vector is divergence-free and curl-free in a domain, and its normal 
component (or tangential components) on the boundary is zero, then this vector is zero. 
This paper emphasizes that one must apply the div and curl operators together with the 
boundary condition (either the normal component or the tangential components be zero) 
to a vector differential equation to derive an equivalent higher-order version of system. 
We believe that many controversies over the permissibility of the boundary conditions for 
derived equations can be resolved via this div-curl method. 

The paper is organized as follows. In Section 1, we use an elementary method (in- 
tegration by parts) to prove Friedrichs’ inequality related to the div and curl operators. 
In Section 2, we show that the three-dimensional div-curl system is properly determined 
and elliptic by introducing a dummy variable, and that the least-squares method for the 
div-curl system corresponds to solving three independent Poisson equations of three veloc- 
ity components with three coupled boundary conditions. In Section 2.4, we introduce the 
div-curl method to change the low-order partial differential equations into an equivalent 
higher-order form. In Section 3, we study the div-curl system with a different boundary 
condition. In Section 4, we use the results obtained in Section 1-3 to prove the theorem 
about the orthogonal decomposition of vectors, and use it to establish Friedrich’s second 
inequality. In Section 5, we show that the Navier-Stokes problem written in the first-order 
velocity-pressure-vorticity formulation is properly determined and elliptic when the con- 
straint condition (the divergence of the vorticity should be zero) is supplied. In Section 6, 
we show that the number of permissible boundary conditions is four for three dimensional 
problems and list all possible combinations of boundary conditions for the Stokes prob- 
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lem. In Section 7 we prove that under four different nonstandard boundary conditions the 
Stokes problem has a unique solution. In Section 8, we derive the second-order velocity- 
pressure-vorticity formulation and the corresponding boundary conditions for the Stokes 
problem via the least-squares method. In Section 9, we derive the conventional second- 
order velocity-vorticity formulation by using the div-curl method. Concluding remarks 
are given in Section 10. In Appendix we list all useful equalities on vector operations, 
briefly derive some important Green’s formulae, and give two Poincare inequalities for 
scalar functions. 


1. Basic Inequalities 


1.1 Notations 


Let f i C 3ft 3 be an open bounded domain with a boundary T, x = z ) a point in 

fl, n be a unit outward normal vector on the boundary, t be a tangential vector to T, and T\ 
and T 2 be two orthogonal tangential vectors to T. In order to pay attention to basic ideas 
and maintain simplicity we often further assume that the domain S7 is simply connected 
and the boundary T is smooth enough, although in many cases these restrictions are not 
necessary, denotes the space of square-integrable functions defined on f l equipped 

with the inner product 


(«>») = / 
J fi 


uvdft u y v £E L 2 (Q) 


and the norm 

IMIo,n = (“>“) ueL 2 {tt). 

denotes the Sobolev space of functions with square-integrable derivatives of order 
up to r. | • | r> n and || • || ri n denote the usual seminorm and norm for H r (Q,), respectively. 
For vector-valued functions u with m components, we have the product spaces 


£ 2 (ft) m , H r (Sl) m 


with the inner product 


(u,v) = / 

J ci 

and the corresponding norm 


u • vdfl u,v € L 2 (fi) T 


Mlo,n — ll u illo,n7 IMIr.n — 51 

j = i j =1 
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Further we define 


< u,v >r= j uvds u E H 1 / 2 ( r), v E H 1 / 2 (r) 


and the corresponding norm 

IMI?/ 2( r =< u ’ u >r • 

When there is no chance for confusion, we will often omit the measure Q or T from the 
inner product and norm designation. 

Througout the paper C denotes a positive constant dependent on ft with possibly 
different values in each appearance. 


Lemma 1.1, Let ft be a bounded open subset of 3ft 3 with a piecewise C 1 boundary I\ 
Every function u of fT^ft) 3 with n x u = 0 on T satisfies 

1*1? + Air + = II v - “Iti + II v x (i.i) 

Jr ^2 

where i?i and R 2 denote the principal radii of curvature for I\ 

Proof. Using Green’s formulae (B.2), (B.3) and (B.7), and Equality (A. 4), we have 

||V.«||J + ||VX5||J 


= (V ■ u, V ■ u) + (V x u, Vxu) 

= ( u , — V(V • u))+ < V • u, n -u > +(u, V(V • u) — Au)+ < V x u, n xu > 

= (ViZ, Vu)+ < V ■ ti, n • tZ > + < V x u, n x u > — < — , u > . (1.2) 

On 

Obviously 

(Vu, Vu) = \u\\. 

Now we turn to the boundary integral terms, n x u = 0 on T implies that 

u = Un on r, 


where U is a scalar function which depends on the location on the boundary surface. By 
using (A.l), the boundary integral terms in (1.2) can be written as follows: 


= J {U{UV -n + Ti-VU)-UVU -n}ds 

-L 


U 2 V -fids. 


It can be verified that on a curved surface 


_ _ 1 1 

V • 71 = — + — = 2 7, 

-ill xi2 


where 7 is the mean curvature. 

□ 


Lemma 1.2. Let be a bounded open subset of 9? 3 with a piecewise C 1 boundary T. 
Every function u of H 1 (J7) 3 with n • u = 0 on T satisfies 

|s|i + j \u 2 is = II V -u||J+ II V x s|i;, (1-3) 

where #1 < R < iZ 2 , in which Ri and R 2 denote the principal radii of curvature for I\ 
Proof. In this case we still have (1.2). Since n • u = 0, we may assume that 

u = Uf on r. 


By virtue of the triple scalar product 

(V x u) • (n x u) = n • (u x (V x u)) 
and using (A.5), the boundary integral terms in (1.2) can be written as 

1 du 2 . 


f r „ _ N 1 du\. 

-)J, 

= n.(V(^ 2 )- (»V)E) - 

= J {n 1 ( — (uV)u)}<is 

= J {U 2 n ■ (-(rV)r)}ds 
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= ■(]?)}*>■ 

So the Lemma is proved. Here R is the radius of curvature of the boundary surface in the 
direction of r . 

U 


Since the curvature 1/R is always positive when the boundary surface is convex, we 
derive immediately the following theorem: 

Theorem 1. Let ft be a bounded and convex open subset of & 2 3 , and (Ei ,1^) be the pieces 
of the piecewise C 1 boundary surface I\ Either T i or T 2 may be empty or with strictly 
positive measure. Every function u of H l (fl) z with n • u = 0 on and n x u = 0 on 1^2 

CO tl cfipc • 

l«l?<l|V-5||i + ||Vx5||g. (1.4) 

Theorem 2. Suppose that the simply connected Cl and T satisfy the same assumption in 
Theorem 1. If u 6 /^(O) 3 satisfies V • u = 0, V x u = 0, n • ulp! — 0 and n x tZ|r 2 = 0, 
then u = 0. 

Proof. From Theorem 1, we have 

M? < o, 

that is, u must be a constant vector in Q. From the boundary conditions we know that 
this constant vector must be zero. 

□ 

Theorem 3 (Friedrichs’ First Inequality). Suppose that the simply connected Cl and T 
satisfy the same assumption in Theorem 1. Every function u of with n • u = 0 on 

Ti and n x u = 0 on satisfies: 

M?<C'(||V.«||j + ||Vx5||S), (1.5) 

where the constant C > 0 depends only on Cl. 

The proof of Theorem 3 can be based on the use of contradiction arguments together 
with Theorem 1 and 2, see e.g., Saranen[21] and the references therein. In the two- 
dimensional case, a direct proof is available[22,23]. 


2. The Div-Curl System (Case 1) 

2.1 Determinedness and ellipticity 
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Let us first consider the following three-dimensional div-curl system: 


V ■ u = 0 in Cl y 

(2.1a) 

V x u = uj in fl, 

(2.16) 

n • u = 0 on I\ 

(2.1a) 


where the domain ft is bounded and convex with a piecewise C 1 boundary T , the vector 
function 6 L 2 (ft) 3 is given and satisfies the compatibility conditions: 

V -W = 0 in ft, (2.2a) 


n • a >ds — 0. 


(2.26) 


The first-order system (2.1) is fundamental for incompressible viscous flow problems in 
which u represents the velocity, and a? the vorticity. The system (2.1) also governs, for 
example, static magnetic problems in which u is the magnetic field, and W the electric 
current density (assume that the permeability is one). 

At first glance, System (2.1) seems overdetermined, since there are four equations and 
three unknowns (i.e. three components of the velocity vector). Surely, for this reason solv- 
ing (2.1) is not easy by conventional finite difference or finite element methods. However, 
after careful investigation by applying the same trick as used in [24], we shall find that 
System (2.1) is properly determined and elliptic. By introducing a dummy variable <f>, 
System (2.1) can be written as 

V -u = 0 in ft, (2.3a) 


V0 + V x u = in ft, 
u • n = 0 on r, 

<f> = 0 on r. 


(2.36) 

(2.3c) 

(2.3d) 


Substituting (2.3b) into (2.2a) and taking into account that V • V x u = 0, we obtain 
A<f> = 0. Since <f> = 0 on T, we must have <f> = 0. Therefore System (2.3) is indeed 
equivalent to System (2.1). In Cartesian coordinates System (2.3) is given as 

du dv dw _ 

dx + dy dz ’ 

dd> dw dv 

dx + dy dz 

d* + du_dw = (2.4) 

dy dz dx 
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d<f> dv du 

dz ^ dx dy z 

Now we write System (2.4) in the standard matrix form: 


in which 


A! = 


Ax 

du 

dx 

+ a 2 

du 

dy 

A & U 

+ As 8J 

+ Au 

= f, 


/I 

0 

0 

0\ 


0 

1 0 

°\ 

1 0 

0 

0 

1 1 

, a 2 = 

0 

0 1 

0 1 

° 

0 

-1 

0 i 

0 

0 0 

1 

Vo 

1 

0 

0 J 


\-l 

0 0 

0/ 


A 3 = 


/° 

0 

1 

°\ 

/° 

0 

0 

°\ 

0 

-1 

0 

0 I . 

0 

0 

0 

0 1 

1 

0 

0 

0 ’ A_ 

0 

0 

0 

0 ’ 

Vo 

0 

0 

1/ 

Vo 

0 

0 

0/ 


f = 


/° 






u = 



(2.5) 


The characteristic polynomial associated with System (2.4) 


dei(Ai£ A-2?7 H“ -^sC) = 


* 

0 

C 

V-77 


*7 

-< 

0 


C OX 

\ i I = ^ + c 2 ) 2 # 0 


for all nonzero triplet (£,77,C)> System (2.4) and thus System (2.1) is indeed elliptic and 
properly determined. 

The fist-order elliptic system (2.3) has four equations and four unknowns, so two 
boundary conditions are needed to make System (2.3) well-posed. Here <j> = 0 serves as 
one boundary condition, and n • u = 0 as another one. 


2.2 The least-squares method 

Now let us consider the least-squares method for solving the div-curl equations (2.1), 
see also a theoretical analysis by Fix and Rose[25] aind an application by Hafez and Soli- 
man[26]. We construct the following quadratic functional: 

I : H -=± ft, 
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J(tt) = ||V.tt||2 + ||Vx«-5;||2, 

where H = {u £ fl’ 1 (fi) 3 |n • u — 0 on T}. We note that the introduction of a dummy 
variable <f> in § 2.1 is only for the verification of the determinedness, and is not required 
in the least-squares functional I. Taking the variation of I with respect to u, and letting 
5u = v and 81 — 0, we obtain the least-squares weak formulation: Find u£H such that 

B(u,v) = L(v) Vv £ H, (2-6) 


where 

B(u,v) = (V • u, V *r) + (V x u, V x v ), 

L(v) = (u>, V x v). 

Obviously B(u,v) is symmetric. Due to Theorem 1 we immediately see that 

|u|i < B(u,u) = L(u) < |u|i||u;||o, 

and 

|u|i < H|o. ( 2 - 7 ) 

If, in addition, f 1 is simply connected, due to Theorem 3 we also see that 

^Mli < B(u,u) - L(u ) < ||u||iN|o. 

By virtue of the Lax-Milgram theorem, in fact we have proved the following Lemma. 

Lemma 2.1 If ffc is bounded, convex and simply connected, then the solution of (2.1) 
uniquely exists and satisfies: 

Nil < C|M|o. (2.8) 

2.3 The Euler-Lagrange equation 

In order to understand how the least-squares method works, we derive the Euler- 
Lagrange equations associated with the least-squares weak formulation (2.6) which can be 
rewritten as: Find u £ H such that 


C V • U, V • v) + (V x u - U>, V X v) = 0 Vv £ H. (2.9) 

Suppose that u and uj are sufficiently smooth. By using Green’s formulae (B.3) and (B.5), 
Equation (2.9) can be written as 

(— V(V • u),v)+ < V -u,fi -v > + 
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( 2 . 10 ) 


(V x (V Xu — u>),v)— < n x (V Xu — w), v >= 0. 

Taking into account (A. 4) and that v satisfies n -v = 0, from (2.10) we obtain 

(— Au — V x u,v)— < n x (V xu — w),v >— 0 (2.11) 

for all admissible v 6 H, hence we have the Euler-Lagrange equation and boundary con- 
ditions: 

— Au = V x u> in ft, (2.12a) 

n x ( V x u — w) = 0 on T, (2.126) 

n • u = 0 on r. (2.12c) 

We remark that included in the boundary integral term in (2.11) is a triple vectqr 
product of n, (V x u — uJ) and v. Since we already know that v is orthogonal to the normal 
n, to make the triple vector product be zero requires only that (V x u-u) is parallel to 
to on r, which is represented algebraically by (2.12b). It is not necessary to require that 
(V x a - w) = 0 on T. 

The first-order div-curl system (2.1), the least-squares weak formulation (2.9), the 
second-order Poisson equations (2.12), and the Galerkin weak form (2.11) are all equiva- 
lent. Now it turns out that the least-squares method (2.9) for the div-curl equations (2.1) 
corresponds to using the Galerkin method (2.11) to solve three independent Poisson equa- 
tions (2.12a) with three coupled boundary conditions (2.12b) and (2.12c). Here (2.12b) 
serves as two natural boundary conditions, and (2.12c) as an essential boundary condition. 
The attractions of using (2.12) are obvious. One avoids dealing with the incompressibility 
constraint (2.1a); instead, one deals with Poisson equations that everyone would rather 
solve. 

It is well known that the Galerkin method is a perfect method for Poisson equations. 
Here “perfect” means that the corresponding finite element method has an optimal rate 
of convergence and leads to a symmetric positive definite matrix. This fact explains why 
the least-squares method is a perfect method for the first-order system (2.1). 

Unfortunately, this perfect least-squares method is often misunderstood. Some people 
think that “the Euler-Lagrange equation derived from the least-squares weak formulation 
is some derivative of the original equation, hence this introduces the possibility of spurious 
solutions if incorrect boundary conditions are used.” However, our derivation clearly shows 
that the least-squares method does not require any additional boundary condition. Only if 
someone would like to use, for example, the finite difference method to solve the associated 
Euler-Lagrange equation (2.12a), additional natural boundary conditions (2.12b) (which 
are simply taken from the original first-order system) are needed in order to obtain the 
solution. 

2.4 The div-curl method 
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The derivation in §2.3 shows that the least-squares method converts the difficult first- 
order div-curl system into an easy second-order system and reveals that each component of 
the velocity satisfies a Poisson equation. Here we show how to derive a high-order version 
of the differential equations without any risk of generating spurious solutions by another 
systematic way. By virtue of Theorem 2, System (2.1) is equivalent to 


V • u = 0 in ft, 

(2.13a) 

V ■ (V x u — u>) = 0 in ft, 

(2.136) 

V x (V x u — ui) — 0 in ft, 

(2.13c) 

nx(Vxu-u>)=0 on T, 

(2.13d) 

n ■ i 1 = 0 on F. 

(2.13e) 

Due to the compatibility constraint (2.2a) and V ■ V x u = 0, 
satisfied. After simplification System (2.13) becomes 

Equation (2.13b) is always 

—Au = Vxw in ft, 

(2.14a) 

n x ( V x u — u;) = 0 on T, 

(2.146) 

n • u = 0 on r, 

(2.14c) 


which is the same as Equation (2.12) obtained by the least-squares method. Here we remark 
that by this div-curl method itself it is not clear whether the divergence-free equation of 
the velocity, i.e. Equation (2.13a), can be eliminated or not. It is only through the least- 
squares method that we are able to make sure that the satisfaction of Equation (2.14a) and 
boundary condition (2.14b) and (2.14c) can guarantee that the solution of u is divergence- 
free. 

Yet from Theorem 2 we may choose another boundary condition to replace (2.14b), 
so that we solve the following system: 


V • u = 0 in Ct y 

(2.15a) 

-Au = V x u> in ft, 

(2.156) 

(V xu-w) = 0 on r, 

(2.15c) 

n • u = 0 on r. 

(2.15d) 
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In this case the divergence-free equation (2.15a) cannot be eliminated. 


3. The Div-Curl System (Case 2) 


Let us consider the following three-dimensional div-curl equations with another im- 
portant boundary condition: 

V-u = p in Cl, (3.1a) 

Vxu = 0 in Cl, (3.16) 

n x u = 0 on T, (3.1c) 

where p is given. The first-order system (3.1) governs, for example, static electric problems 
in which u is the electric field, and p the density of the electric charge. System (2.1) has only 
one algebraic boundary condition, while System (3.1) includes two independent algebraic 
boundary conditions. We shall explain why these two boundary conditions are correct. 

At first we show that System (3.1) is also properly determined and elliptic. By intro- 
ducing a dummy variable <f>, System (3.1) can be written as 


V • u = p in Cl, (3.3a) 

V^+Vx« = 0 in Cl, (3.36) 

n x u = 0 on T. (3.3c) 

Notice that in this case we don’t impose any boundary condition for the dummy variable 
<f>. By virtue of Theorem 2, (3.3b) is equivalent to the following equations and boundary 
condition: 

V ■ (V<^ + V x u) = 0 in 11, (3.4a) 

V x (V<f> + V x a) = 0 in fl, (3.46) 

i'(V^|Vxu) = 0 on T. (3.4c) 

Since n x u = 0 on T, that is, u is parallel with n on T, or u is perpendicular to T, we 
necessarily have 

n ■ (V x«) = 0 on T. (3-5) 

The proof of (3.5) is straightforward. Assume the contrary, say, n • (V x u) > 0 at a point 
P on T, then in a neighbourhood dr of P we have 


n • (V x u) > e > 0, 
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in which e is a small positive constant. From the Stokes theorem we have a contradiction: 




• nds > 0, 


where c is the boundary contour of dT. 

Taking into account (3.5), (3.4a) and (3.4c) lead to 


A<f> = 0 in 0, 

(3.6) 

^ = ° on T. 
dn 

(3.7) 


That is, <f> is a constant, or V<f> = 0 in fl. Therefore System (3.3) with four equations 
and four unknowns is equivalent to System (3.1). System (3.3) is elliptic and properly 
determined, so is System (3.1). The boundary condition (3.1c) requires that two tangential 
components of u must be zero on the surface T. In this case no boundary condition on 
the du mm y variable <j> should be specified, so altogether there are only two boundary 
conditions which are consistent with a 2 x 2 first-order elliptic system. 


4. Orthogonal Decomposition of Vectors 
Theorem 4. Every vector u G 5 1 (fl) 3 has the orthogonal decomposition: 

u — Vg + V x ip, (4-1) 

where q G jff 2 (0)/Jf and ip G H 2 (Cl) 3 . 

Proof, we note that for the purpose of investigation in this paper, there is no need to 
find q and ip. This theorem shall be proved by directly seeking for Vg and V x ip. By 
virtue of Theorem 2, Equation (4.1) is equivalent to the following equations and boundary 


condition: 

V • (Vg + V x if — S) = 0 in fi, (4.2a) 

V x C Vg + V x - u) = 0 in ft, (4.26) 

to • ( Vg + V x rp — u) = 0 on T. (4.2c) 

Taking into account V • V x ij> = 0 and V x Vg = 0, System (4.2) can be written as follows: 

Ag = V ■ u in fl, (4.3a) 

V x (V x ^) = V x u in ft, (4.36) 
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n ' (Vq + V x \j>) = n • u on I\ (4.3c) 

We may solve the following Neuman problem of Poisson equation to obtain q: 

A q = V • u in ft, (4.4a) 

n • Vq = n • u on T, (4.46) 

where the boundary condition (4.4b) is additionally supplied. Although q is not unique, 
i.e., an arbitrary constant can be added into q, Vq is uniquely determined. 

Now t/> should satisfy 

V • (V x ip) = 0 in ft, (4.5a) 

V x (V x VO = V x u m ft, (4.56) 

n • (V x V) = 0 on r. (4.5c) 

V x V in System (4.5) may be considered as an unknown vector and can be uniquely 
determined by the least-squares method described in Section 2. Therefore the validation 
of the decomposition (4.1) is proved. 

Using (B.4) and (4.5c) we have 

(Vq, V X V 1 ) =< n • (V X ^),q >= 0, 
that is, Vq and V x ^ are orthogonal. 

□ 

Since q is the solution of the Neuman problem of Poisson equation (4.4), we have the 
following regularity result [27,28]: 

|qk < C{||V • u||o + ||n • tt||i/ 2 ,r}- (4-6) 

Since V x ip is the least-squares solution of (4.5), by virtue of (2.7) we have 

|V x V>|i < ||V x u|| 0 . (4.7) 

Hence by using (4.1), (4.6) and (4.7) we obtain 

INIi < l|Vg||i + ||V x ^||o + |V x ^ 

< ||Vq||i + ||u]| 0 + ||Vq|| 0 + |V x 
< ||u||o + i ql 2 + |V x ifh 
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< C’{||tt||o,n + || V • u||o,n + ||« • «lli/ 2 ,r + ll^ 7 x «llo,n}- 

In fact we have proved the following Theorem. 

Theorem 5 (Friedrich’s Second Inequality). Let Cl be a bounded and convex open region 
of with a piecewise C 1 boundary T. Every u € satisfies 

llujli.n < C'{|| : u||o,n + ||V • u||o,n + II V x u||o,n + ||« • «||i/ 2 ,r}* (4-8) 


5. The Velocity-Pressure- Vorticity Formulation 
of the Navier-Stokes Equations 


Usually the Navier-Stokes equations governing the steady-state incompressible viscous 
flows are written in the following velocity-pressure formulation 


u • Vu + Vp — — Au = / 
Re 


V • u = 0 


in Q, 


171 O* 


(5.1a) 

(5.16) 


We assume that the domain Cl is bounded, convex and simply connected with a piecewise 
C 1 boundary T. All variables in (5.1) are nondimensionalized, u = ( u,v,w ) denotes the 
velocity, p the pressure, / = (/*, / 9 , fz) € L 2 (Clf the body force, and Re the Reynolds 
number, defined as 

UL 

Re = , 


v 


where L is a reference length, U a reference velocity and u the kinematic viscosity. 

All mathematical analyses of the existence of the solution of (5.1) are conducted by 
using the Galerkin method mainly under the standard velocity boundary condition 

u = up 071 r, (5.1c) 


where «r denotes a given function defined on the boundary T . 

It is well known that the Galerkin mixed method leads to a saddle-point problem, 
thus the sophisticated LBB condition is invoked to guarantee the existence of the solution. 
It is notoriously difficult to verify and satisfy the LBB condition. From a numerical point 
of view, the most difficult problem associated with the Galerkin mixed method is that the 
resulting discretized algebraic equations are nonsymmetric and nonpositive-definite which 
are hard to deal with for large problems. All these difficulties motivated us to apply the 
least-squares method. 
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Since the momentum equation (5.1a) involves the second-order derivatives of veloc- 
ity, the application of the least-squares method to (5.1) requires the use of impractical 
continuously differentiable functions. In order to avoid this trouble, one has to consider 
the governing equations of incompressible flow in the form of a first-order system. The 
velocity-pressure-stress formulation is one of the choices. However, this formulation has 
too many unknowns. Moreover, as pointed out in [15], the three-dimensional velocity- 
pressure-stress formulation has nine independent unknowns and nine independent equa- 
tions, and thus cannot be elliptic in the ordinary sense. Consequently, the selection of 
proper boundary conditions becomes a difficult task. Instead, we introduce the vorticity 
u; = , u>j, , u } z ) = V x u as an additional independent unknown vector, and rewrite the 

incompressible Navier-Stokes equations in the following first-order quasi -linear velocity- 
pressure-vorticity formulation: 


u • Vu -I- Vp 4- ■— V x i») = / in ft, (5.2a) 

Re 

V • Ta = 0 in ft, (5.26) 

a; — V x u = 0 in fl, (5.2c) 

V • u = 0 in ft. (5. 2d) 


Here we have included the compatibility constraint condition (5.2b), i.e., the divergence 
of the vorticity vector equals zero, to make System (5.2) elliptic in the ordinary sense. 
The determinedness and ellipticity of (5.2) have been proved in [15]. For completeness we 
briefly repeat the proof in the following. 

At first glance, one may think that System (5.2) is overdetermined, since there are 
seven known variables, i.e., three velocity components u,v,w, one pressure p and three 
vorticity components iv x ,u> y ,u> z , and eight equations. We shall show that System (5.2) is 
really properly determined and elliptic by using the same technique as discussed in Section 
2. As explained in [15], the nonlinear convective terms and the Reynolds number have no 
effect on the classification, so we may just consider the following first-order system of the 
Stokes problem: 

Vp + V x u> = / in ft, (5.3a) 

V • a; = 0 in Cl, (5.3 b) 

—uj + V<f> + V x u = 0 in fl, (5.3c) 

V ■ u — 0 in fl. (5. 3d) 

Here we have already introduced a dummy variable <f> in (5.3c), which satisfies the boundary 
condition <f> — 0 on T. Substituting (5.3c) into (5.3b) yields A 4> — 0, thus <£ = 0 in Cl. 
That is, the introduction of (f> does not change anything. However, now there are eight 


18 



unknowns and eight equations in (5.3), and System (5.3) and hence System (5.2) is indeed 
properly determined. 

We note that in some cases the specification of the boundary condition for the dummy 
variable <fi is unnecessary, and V</> == 0 can also be guaranteed, see Section 6. 

Now let us classify System (5.3). In Cartesian co-ordinates, System (5.3) is given as 


dp 

i 

dui z 

dW y 

dx 


dy 

dz 

dp 

i 

dw x 

dui z 

dy 

t 

dz 

dx 

dp 

i 

dWy 

dui x 

dz 

-+■ 

dx 

dy 

dui x 


duty 

dui z 

dx 

+ 

dy 

+ dz 


= /*, 

— fyi 
~ fz, 
= 0 , 


d<f> dw dv 

d<f> du dw 

d<f> dv du 

~ u ‘ + Tz + Tx ~ ai = ’ 

du dv dw _ 

dx dy + dz 

We may write (5.4) in the standard matrix form: 


. d\l C?U 0U . - 

Al & + 2 % + ’ 


(5.4) 


(5.5) 


in which 


/0 0 0 000 0 1 

00 0 000 -1 0 

00 0 001 00 
00 0 010 00 
00010000 
00 -1 000 0 0 

01 000000 
VlO 0 000 0 0 


/ 0 000 0 01 0 \ 

0 000 0 001 
0 000 -1 000 
0 000 0 100 
0 010 0 000 ’ 
0 001 0 000 
-1 000 0 000 , 
Vo 100 0 00 0/ 
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/O 0 000-1 0 0\ 

00001 000 
00000001 
_0 0 0000 10 
3 “ 0 -1 000 0 00 ’ 
1 0000000 
00010000 
Vo 0 1 0 0 0 0 0 / 

/f*\ 

fv 

u 

f= 0 

0 ’ 

0 

0 

Vo/ 


/O 0 0 0 0 0 0 0 \ 

0000 0 0 0 0 

00000 0 00 

0 0000 0 00 

0000 -1 0 0 0’ 

0000 0 -1 0 0 

0000 0 0 -1 0 

Vo 0 0 0 0 0 00/ 



The characteristic polynomial associated with System ( 5 . 5 ) 

/ 0 0 0 0 0 — £ 77 £ \ 

0 0 0 0 < 0 — ^ 77 

0 0 0 0 —77 ^ 0C 

det(A 1 Z + A 2 T ) +A 3 () = det jj _°^ ® ® g q q = ifW + C 2 ) 4 ^ 0 

C 0 t] 0 0 00 

-77 £ 0C0 0 00 

V £ V COO 0 0 0/ 

for all nonzero triplet (£, 77, C), System ( 5 . 5 ) or ( 5 . 4 ) and thus System (5.3) is indeed elliptic. 
The ellipticity can be easily understood in the following way. Equations (5.3a) and (5.3b) 
constitute a div-curl system of the vorticity, and Equations (5.3c) and (5.3d) constitute 
a div-curl system of the velocity. These two div-curl systems are coupled through the 
vorticity in Equation ( 5 . 3 c) and form the Stokes equations. In other words, the principal 
part of the Stokes operator consists of two identical elliptic div-curl operators. 

The classification in this paper is based on an ordinary method, so Equation (5.2b) 
is needed to guarantee ellipticity. System (5.2) is also elliptic in the Agmon-Douglis- 
Nirenberg(ADN) sense[ 29 ]. In fact, System (5.2) without (5.2b) is also elliptic in the ADN 
sense. In the ADN theory the non-principal part of the operator must be considered. The 
corresponding least-squares method can be developed based on minimizing a weighted 
(mesh-dependent) L 2 norm of the residuals following the idea proposed in [ 30 ] for general 
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elliptic systems and the Stokes equations in the velocity-pressure formulation, and in [31] 
for the Stokes equations in the velocity-pressure-vorticity formulation. 


6. Boundary Conditions for the Navier-Stokes Equations 

The boundary conditions for the Navier-Stokes equations have been extensively dis- 
cussed in the literature. The most complete and thorough examination of this topic may 
be found in Gresho’s state-of-the-art reports[2,32,33] and the references cited therein. For 
boundary conditions other than the specification of the velocity one may also consult 
Gunzburger[9], Girault[34], Pironneau[10] and Verfurth[35|. Based on different choices for 
the formulation of the viscous term in the velocity-pressure formulation and the Galerkin 
method, Gunzburger correctly gives many possible combinations of nonstandard boundary 
conditions. Here we should mention that for nonstandard boundary conditions there were 
very few rigorous analyses available. In this section we shall list all possible combinations 
of boundary conditions based on the first-order velocity-pressure-vorticity formulation. 

Before investigation of boundary conditions we first discuss the equivalence between 
the velocity-pressure formulation (5.1) and the velocity-pressure-vorticity formulation 
(5.2). System (5.2) is obtained from (5.1) by introducing the vorticity To which is some 
kind of derivative of the velocity u and thus reducing the order of differential operator. Re- 
ducing the order of differential equations in this way does not generate spurious solutions. 
It means that System (5.2) is equivalent to System (5.1), that is, the solution of (5.2) is the 
solution of (5.1). Conversely, we may think that System (5.1) is deduced from System (5.2) 
when the definition of the vorticity (5.2c) is substituted into (5.2a) and (5.2b). This type 
of substitution and combination does not generate spurious solutions either. Therefore, 
the velocity-pressure formulation (5.1) and the velocity-pressure-vorticity (5.2) formulation 
are mutually equivalent. This equivalence implies that the permissible boundary condi- 
tions for System (5.2) must be the permissible boundary conditions for System (5.1). The 
reverse is also true. Consequently, it is not possible that some boundary conditions that 
have been shown to be legitimate for System (5.1) might not be so for System (5.2), and 
vice versa. 

As pointed out in the previous section that the nonlinear convective term has no effect 
on the classification of the Navier-Stokes equations, hence the boundary conditions for the 
Stokes equations are valid for the Navier-Stokes equations. So we need only to analyse the 
Stokes problem (5.3). Since the system (5.3) is of first-order, the boundary conditions do 
not involve the derivatives of unknowns. In other words, there are only essential boundary 
conditions for the solution of first-order partial differential equations. This fact precludes 
the mathematical legitimacy of taking the derivative of pressure as a boundary condition. 

For convenience we rewrite the Stokes problem (5.3) as the following two coupled 
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Table 1 . Combination of Boundary Conditions 


Boundary Conditions 

3D 

2D 

Remarks 

(1) Normal velocity 

n • u 

n • u 


Tangential vorticity 

n XU 

w 


(2) Pressure 

p 

p 


Normal vorticity 

n * U 



Normal velocity 

n • u 

n • u 


(3) Pressure 

p 

mm 


Normal vorticity 

n • U 

mu 


Tangential velocity 

n Xu 



(4) Tangential velocity 

n Xu 



Tangential vorticity 

n x U 

■H 


(5) Velocity 

■ 






(6) Pressure 

p 

p 


Normal vorticity 

n • U 


on a 

Tangential vorticity 

n X U 

to 

part of V 
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div-curl systems: 


Vp + V x w = f in ft. 


(6.1a) 

(6.16) 


V • (jj — 0 zn 


and 

V<j> + V xu=u } in Cl, (6.2 a) 

V • u = 0 in ft. (6.26) 

If the vorticity v in (6.2a) is known, the structure of these two div-curl systems would be 
identical. In Section 2 we have thoroughly investigated this type of div-curl system. In 
order to solve (6.1) to obtain p and w, on the boundary we should specify 

(p and n ■ u?) or n x w. 

When v is obtained, to obtain <f> and u we may solve (6.2) with the boundary condition of 

[<f> and n • u) or n x u. 

From the above consideration we can immediately list four permissible combinations of 
boundary conditions (l)-(4) in Table 1 for the Stokes problem. In Table 1 for three- 
dimensional problems we don’t explicitly include the boundary condition (j> = 0 on T. If 
we understand that <f> — 0 always comes with the condition of n-u and count this condition, 
then the total number of boundary conditions is four for 3D problems. Since there are 
eight equations and eight unknowns in the first-order elliptic system (5.3) or (6.1) and 
(6.2), we need four boundary conditions on each fixed boundary. In (5) and (6) of Table 1 
we list other possible choices which are not from the above consideration but still satisfy 
the requirement of four conditions on the boundary. For two-dimensional problems in the 
first-order velocity-pressure-vorticity formulation there are four unknowns, i.e., u,v,p,uj, 
and four equations, i.e., two momentum equations, one definition of the vorticity and the 
incompressibility; and no dummy variable is involved, see [12]. Therefore, two boundary 
conditions are needed for two-dimensional problems. In Table 1 we also list corresponding 
boundary conditions for two-dimensional problems. 


7. Permisibility of the boundary conditions 

In the this section we shall rigoriously prove the well-posedness of the boundary con- 
ditions (l)-(4) in Table 1. The boundary conditions (l)-(5) in Table 1 can be used on the 
entire boundary or on a part of boundary T. For simplicity we consider only one kind of 
homogeneous boundary conditions on the entire boundary. The results can be extended 
to mixed and inhomogeneous cases without difficulty. 
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Given the elliptic differential operator, the question of well-posedness reduces to ver- 
ification of the permissibility of the boundary conditions. For general elliptic systems one 
may use the ADN theory to accomplish this task. For first-order systems in the plane, one 
may also try to answer this question by applying the theory developed in [36]. However, 
both techniques invoke the modern theory of elliptic partial differential equations and are 
quite difficult to understand by engineers. In this paper we try to identify the permissible 
boundary conditions by using an elementary treatment. The mathematical tools used are 
the least-squares method, Green’s formulae (integration by parts) and Friedrich’s inequal- 
ities established for div and curl operators. An elliptic system with supplied boundary 
conditions is considered to be well-posed, if one can prove that the corresponding least- 
squares method leads to a coersive bilinear form. In the following, we do this case by 
case. 


7.1 u n — 0, w T i — 0, u > r 2 = 0 (n • u = 0, n x w = 0) on T 


These boundary conditions may be used for the symmetric plane. The inhomogeneous 
version may be used for the uniform inflow boundary in which the normal components of the 
velocity and two tangential components of the vorticity are prescribed. These conditions 
correspond to those in the velocity-pressure formulation, i.e., the normal velocity and 
the tangential stresses are given. For example, let us consider a piece of boundary with 
n = (1,0,0). We have that 

u = 0, 


oj 


y 


U) z 


du dw 

!h~!h =0, 

dv du 


From u = 0 on this boundary we deduce that 


Therefore, 


This implies that 


du 

dl~ 0i 

du n 


dw 

dx 

dv 

dx 


= 0, 


= 0. 


du dw 

Ixz = 7) I" = 
dz dx 
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dv du 

lxy = dx + ~dy = °* 

That is, the tangential strains and thus the tangential stresses are zero. 

In order to guarantee the uniqueness of the solution of pressure p, we require that the 
pressure has zero mean over Cl: 

f pdx = 0. (7-1) 

J n 

The least-squares method minimizes the following functional 

J{ U) = || Vp + V x w - /Ho + II V • w||J + ||w - V x u\\l + || V ■ u\\l, (7.2) 

where U = (u,p,u) E H = and U satisfies the corresponding homogeneous 

boundary conditions on T. Furthermore, / E Following standard arguments of 

variational calculus, we deduce that the least-squares weak solution U necessarily satisfies 

B(U,V) = L(V) W = (u, 9 ,r) EH, (7.3) 

in which 

B{ U,V) = (Vp + V xui,Vg + Vxf) + (V-«,V-f)+ 

(a; — V x u,t — V x v) 4- (V • u, V • u), (7-4) 

i(V) = (/,Vg + Vxf), (7.5) 

and V satisfies the same homogeneous boundary conditions as U. 

Clearly i?(U, V) is symmetric. If we can prove f?(V, V) is coersive, then the existence 
and the uniqueness of the weak solution follow from the Lax-Milgram theorem in a standard 
m ann er. Consequently, the corresponding finite element method has an optimal rate of 
convergence for all unknowns. 

Now we examine the coersivity of i?(V, V). We have 

B{V, V) = || Vq + V x r\\l + ||V • t||q + ||t - V x v\\ 2 0 + || V ■ w||J. (7.6) 

Let us expand the first term in (7.6). Since n x r = 0 on T, using Green’s formula (B.6), 
we have 

( Vg, V x r) =< Vg, n x r >= 0, (7.7) 
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and thus 

II V, + V x r|| l = ||V,||? + ||V x f|| l + 2(V S , V x f ) = ||V,||? + ||V x f||?. (7.8) 

By virtue of (7.8) we have 

B(V,V) = ||Vj||J + II V x f||? + ||V .r||5 + ||f - V x *||> + ||v.5||». (7.9) 

From (7.9) we have 

B(V,V)>||V,||S = |,i;. (7.10) 

Since q satisfies the zero mean constraint (7.1), from the Poincare inequality (C.l) we have 

C\q\\ > h\\l ( 7 . 11 ) 

Combining (7.10) with (7.11) yields 

CB(V,V) > ||,||?. (7.12) 

From (7.9) we also have 

B(V,V)>||Vxf||i + ||V.f||?, (7.13) 

S(V,V)>||f-Vx»||5, (7.14) 

«(V,V)>||V.5J|J. (7.15) 

Since n x r = 0 on F, from Theorem 3 (Friedrichs’ inequality) we have the inequality: 

car v x f ii? + ii v ■ f n?) > nr ii? > Hr n§. (7.i6) 

Combining (7.13) with (7.16) yields 

CB(V,V)> ||r||? (7.17) 

and 


or 


From (7.14) we have 


CB(V,V) > ||f||J 
C(B(V,V))i > ||f|| 0 . 
(B(V,V))i >||-f + Vx»||„. 


(7.18) 

(7.19) 
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Combining (7.18) with (7.19) and using the triangle inequality we have 
C(B(V,V))* > || t || 0 + || -r + Vx v|| 0 > ||V x t?|| 0 , 

that is 

C(B(V,V)) > || V x u||q. (7.20) 

Combining (7.15) with (7.20) leads to 

CB(V,V) > ||V x v\\l + ||V -v\\l (7.21) 

Since n • v = 0 on T, again from Theorem 3 we have the inequality 

(7(|| V x v||q + || V *v||o) > IMIi* (7-22) 

Combining (7.21) with (7.22) yields 

CB(V,V) > ||»||?. (7-23) 

Combining (7.12), (7.17) and (7.23) together we finally obtain that 

CB(V,\) > ||v||i + ||q||i + |Mli- (7-24) 

This shows that 5(V, V) is indeed bounded below in H 1 norm and thus coersive. Con- 
sequently, it is trivial to prove that this problem has a unique solution that satisfies the 
following bound: 

INI. + IIpII. + Ml s cm- 

7.2 p — 0, U n — 0, U3n = 0 (p = 0, n . u = 0, n • To = 0) on T 

The related inhomogeneous case represents, for example, the well developed inflow 
boundary, in which the normal velocity is given, and the normal vorticity and the pressure 
are prescribed be zero. In two-dimensional cases, only u n and p are prescribed. These 
boundary conditions seem difficult to justify by the Galerkin method. The numerical 
results can be found in Bochev and Gunzburger[19,3l]. 

The corresponding least-squares method minimizes the same functional (7.2). The 
proof of the coersivity of B(V, V) follows the steps of the previous case. 

We have 

B(V,V) = ||Vg + V x r||o + ||V • t||q + ||r — V x v\\l + ||' V • v\\ 2 0 . (7.25) 
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We expand the first term in (7.25). Since g = 0 on T, using Green’s formula (B.4), we have 

(Vg, V x r) ==< g, n • (V x f) >= 0, (7.26) 

and thus 

II Vj + V x T||| = ||Vj||| + II V x T||| + 2( Vj, V x r) = ||Vj||J + ||V x T|||. (7.27) 

By virtue of (7.27) we have 

S(V, V) = II V, 115 + I! V x r||| + ||V • Till + ||T - V x «n: + ||V •®|||. (7.28) 

From (7.28) we have 

B(V f V)>||V®||S = |g|;. (7.29) 

Since q = 0 on T, from Poincare inequality (C.2) we have 

C\q\\ > ||?||?. (7.30) 

Combining (7.29) with (7.30) yields 

C-B(V.V) > Hsllf. (7.31) 


From (7.28) we also have 


•B(V,V) > j|V X r||| + ||V- r|||, 

*(V,V)>||T-VxS|||, 

«(V.V)>||V.5|||. 

Since n • r = 0 on from Theorem 3 we have the inequality: 

C(||' V X r||| + || V- t|||) > ||r||| > ||r|||. 
Combining (7.32) with (7.35) we have 

CB(V,V) >||r||? 

and 

CB(V,\) > ||r||| 

or 

C(5(V,V))I > ||r|| 0 . 


(7.32) 

(7.33) 

(7.34) 

(7.35) 

(7.36) 


(7.37) 
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From (7.33) we have 


(7.38) 


(B(V,V))5 > || -T + V xu||o. 

Combining (7.37) with (7.38) and using the triangle inequality lead to 
C(B(V, V)) * > ||r|| 0 + || - r + V x v|| 0 > ||V x t>|| 0 , 

that is 

C(B(V,V))>||Vx*||j|. (7.39) 

Combining (7.34) with (7.39) leads to 

CB(V,V) > ||V x t;||j| + ||V -v\\l. (7.40) 

Since n • v = 0 on T, from Theorem 3 we have the inequality 

(7(|| V x u||o + || V -u||o) > ll®lli- ( 7 - 41 ) 

Combining (7.40) with (7.41) yields 

C5(V,V)>||v||?. (7.42) 

Combining (7.31), (7.36) and (7.42) together we finally obtain 

c7B(v,v)>ii5|i; + u i |i; + iir|i;. ( 7 . 43 ) 

Therefore, the coersivity of 2?(V, V) is proved. 

7.3 p = 0, u T i = 0, u T 2 = 0, u; n = 0 (p = 0,n x u = 0, n • w = 0) on T 

This boundary condition may be used for the well developed exit boundary. Here four 
boundary conditions are prescribed. As mentioned in §3.3 n x u = 0 on T analytically 
implies that n ■ (V x u) = 0 on T. It seems that there are too many boundary conditions. 
In the previous cases we have specified the boundary condition <f> — 0 on T for the dummy 
variable <f> in advance, so only three boundary conditions are needed. In the following we 
show that in the present case no boundary condition is needed for the dummy variable <f>, 
so it is all right to specify four conditions. 

By virtue of Theorem 2, Equation (5.3c) is equivalent to the following equations and 
boundary condition: 

V x (— To + V<f> + V x u) = 0 in ft, (7 .44a) 

V • (—w + V0 + V x u) = 0 in fl, (7.446) 

n • (— + V4> + V x u) = 0 on T. (7.44c) 
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Taking into account (5.3b) and V • V x u = 0, Equation (7.44b) becomes 


A <f> = 0 in SI. (7.45a) 

Taking into account w n = 0 on T and (3.5), the boundary condition (7.44c) becomes 

(7.456) 


d<f> n 

— = 0 on T. 
an 


Equation (7.45a) and (7.45b) imply that <j> is a constant or = 0 in SI. Therefore, four 
conditions in the present case automatically guarantee that the dummy variable <j> can be 
eliminated in Equation (5.3c). 

These boundary conditions correspond to those in the velocity-pressure formulation, 
i.e., the tangential velocity components and the normal stress are prescribed. To show this 
let’s consider, for example, the surface with n = (1,0,0). Since v = w — 0, we have 


dv 

dy 

dw 

dz 


= 0 , 


= 0 . 


Hence from the continuity of velocity we know that 

du 


dx 


= 0 . 


Therefore, 


that is, the normal stress is zero. 


du 

(T x =p + 2v— = 0, 

Ox 


The least-squares method minimizes the same functional (7.2). We shall now prove 
the coersivity of B(V, V). We have 

B(V, V) = II V, + V x r||| + II V ■ Till + ||r - V x ®||| + ||V • 5|||. (7.46) 

Since q = 0 on T, by virtue of (7.26) we have 

II v 9 + V x r||| = ||V,||| + II V x t|||. (7.47) 

Thus we obtain 

B(V, V) = IIVjIII + II V x Till + || V • Till + ||t - V x s||| + II V • t|||. (7.48) 
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Therefore 


S(V,V)>||V«||2 = |«li- 

(7.49) 

Since q — 0 on T, from Poincare inequality (C.2) we have 

C\q\i > Nil- 

(7.50) 

Combining (7.49) with (7.50) yields 

C£(V,V) > llflft. 

(7.51) 

From (7.48) we know that 

B(V,V) > ||V xr||5 + ||V-r||J. 

(7.52) 

Since n ■ r — 0 on T, from Theorem 3 we have the inequality 

C(||Vx Till + ||V- Till) >||rl||. 

(7.53) 

Combining (7.52) with (7.53) yields 

CB(V,V) > ||r||?. 

(7.54) 

From (7.48) we also know that 

B(V,V) > ||x - V x nil?, 

(7.55) 

B(V,V) >||V.»||J. 

(7.56) 

From (7.54) we have 

CB(V,V)i > ||r||„. 

(7.57) 

From (7.55) we know 

£(V, V)5 > || — r + V x v|| 0 . 

(7.58) 

Combining (7.57) with (7.58) and using the triangle inequality lead to 
CB{V, V)s > ||r|| 0 + || - r + V x v|| 0 > ||V x v|| 0 , 

that is 

CU?(V,V) > || V x v||q- 

(7.59) 

Combining (7.59) with (7.56) yields 

CB(V,"V) > || V x t?||o + || V • v||q. 

(7.60) 
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(7-61) 


Since n x v = 0 on I\ from Theorem 3 we have the inequality 

C(I|V x v||o + || V • v||q) > 

Combining (7.60) and (7.61) leads to 

CB(V,V)>||S||;. (7.62) 

Combining (7.51), (7.54) and (7.62) together we finally obtain the coersivity 

CB(V,V) > Ml? + ||®||f + ||r||;. (7.63) 

If we really specify the dummy variable <f> be zero on T in advance, then only three 
boundary conditions are needed, and that cv n = 0 can be imposed in a weak sense. In this 
case, the least-squares method minimizes the following functional: 

J(U) = ||Vp+ V x S-7II5 + ||V -SllS + ||S - V X s||» + ||V -u\\l + \\n-WWiw (7.64) 
We have 

.B(U, V) = (Vp + V x u>, V 5 + V xr) + (V-w,V'r) 

+(u> — V x u,t — V xu) + (V-u,V •»)+ < W ■ n,r • n >, (7.65) 

i(V) = (/,V, + VXT), (7.66) 

and 

B(V, V) = II V, + V x f 115 + II V ■ r||o + ||r - V x »||5 + ||V ■ S||S + ||n • r||f /2ir . (7.67) 

Since q = 0 on T, by virtue of (7.26) we have 

||V ? + V x Till = ||V 4 ||5 + ||V x f||5. (7.68) 

Thus we obtain 

B(V, V) = II V,||| + II V x t|| 5 + ||V • r||5 + ||r - V x 6 ||S + ||V ■ 5|| J + ||B • t||? / 2J .. (7.69) 
Therefore 

B(V,V)>||V,||5=| { |?. (7.70) 

Since 5 = 0 on T, from the Poincare inequality (C. 2 ) we have 

C\q\i > Nil. (7.71) 
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Combining (7.70) with (7.71) yields 

CB(V,V) > HjIIJ. 

(7.72) 

From (7.69) we also know that 

B(V , V) > ||V x r\\l + ||V • r\\l + ||n ■ r\\ 2 1/2X , 

(7.73) 

B(V,V) > ||r — V x v||o, 

(7.74) 

B(V,V) > ||V • w||2- 

(7.75) 

Now let us expand the fourth term in (7.69). Since n x v = 0, using Green s 

formula (B.5) 

we have 

||r - V x v||o = ||r||o + || V x v\\l - 2(r, V x n) = ||t||q + ||V x u||o - 2(V 

x r,n). (7.76) 

Therefore (7.69) becomes 

S(V, V) = \\Vq\\l + ||V X t\\1 + ||V -t\\1 + ||t||o + ||V x v\\* 
— 2(V X r, v) + II V • u||o + ||S • r||j|/ 2tr . 

(7.77) 

Since n x t; = 0 on T, from Theorem 3 we have the inequality: 

C(||V X v||o + || V • v||o) > ||u||i > ||v||o- 

(7.78) 

Multiplying (7.77) by C 2 and taking into account (7.78) we have 

C 2 B(V,V) > C 2 \\r\\l + C 2 ( ||V x v\\l + ||V -v\\l) - 2 C 2 (V x r 


or 

C 2 jB(V,V) > C 2 ||r||o + C||v||o — 2C 2 (V x t,v). 

(7.79) 

From (7.73) we have 

C 3 £(V,V) > C 3 ||V xr||J. 

(7.80) 

Adding (7.79) and (7.80) yields 

(C 3 + C 2 )B(V,V) > C 2 \\r\\l + C\\CV XT- n|| 3 . 

(7.81) 

From (7.81) obviously we have 

CB(V,V) > ||r||J. 

(7.82) 
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Combining (7.73) with (7.82) leads to 


CB(V,V) > ||V x r\\l + IIV • f IIS + ||r||S + ||t ■ n||? /2J . (7.83) 

From Theorem 5 we know that r satisfies 

C{||f||o + IIV -rllo + II V X r||o + ||r -H|| 1/2 j-> > M,. (7.84) 

Combining (7.83) with (7.84) leads to 

CB(V,V) > ||r||f. (7.85) 

From (7.82) we have 

CB(V,V )* > ||r|| 0 . (7.86) 

From (7.74) we know 

CB(V,V)^ > ||-r + Vxn||o. (7.87) 


Combining (7.86) with (7.87) and using the triangle inequality lead to 
CB(V, V)* > ||r|| 0 + || - r + V x v|| 0 > ||V x v|| 0 , 

that is 

CS(V,V)>||Vxr||J. (7.88) 

Combining (7.75) with (7.88) and considering (7.78) yield 

CB{V , V) > ||w||j. (7.89) 

Combining (7.72), (7.85) and (7.89) together we finally obtain the coersivity 

CB(V,V)>|| 4 ||? + |H|; + ||T||f. (7.90) 

7.4 u T i = 0, u T 2 = 0, u; r j — 0, <jJ T 2 = 0 (n x u = 0, n x w = 0) on T 

For the same reason as explained in §7.3, in this case <f> = 0 is guaranteed even no 
boundary condition of (f> is specified. The coersivity of B can be proved by just following 
the steps in §7.1. 

7.5 u n = 0, u T i = 0, u T 2 = 0 (n • tl = 0, n x u = 0) on T 

Obviously this is a standard permissible boundary condition. However, it seems that 
the permissibility cannot be proved by the elementary method presented in this paper. 
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Fortunately, one may rely on the ADN theory to fulfil the task, see Bochev and Gun- 
zburger[31]. 

7.6 p = 0, oj n = 0, o> Tl = 0, u; t2 = 0 (p = 0, n ■ u> = 0, n x « = 0) on T 

Using the boundary conditions p — 0, = 0 one can solve Equation (6.1) to obtain 
p and uj. However, u cannot be uniquely determined by solving (8.3f) with the natural 
boundary conditions in (8.3g) and (8.3f) (see the discussion in the next section). Therefore, 
this combination can only be used on a part of boundary. 


8. Euler-Lagrange Equations Associated with 
the Least-Squares Method for the Stokes Equations 

In order to understand how the least-squares method works, we derive the Euler- 
Lagrange equations associated with the least-squares weak formulation (7.3) for the Stokes 
problems which can be rewritten as 

(Vp + V X a > — /, Vg 4- V X t) + (V *u>, V • r)+ 

(u; - V xu, r - V x u) + ( V • u, V • w) = 0, 

or _ 

(Vp 4- V xu- f, Vq) 

4-( Vp + Vxw — /, Vxt) 

+(V • u?, V • r) 

4-(u> — V X u, t) 

—(uj — V xu, V x u) 

+(V-u, V-v) = 0. (8.1) 

Using Green’s formulae (B.l), (B.3) and (B.5) from (8.1) we have that 

—(V • (Vp 4 - V x a » — /), 9 )+ < n • (Vp 4- V x a ; — /), q > 

4-(V x (Vp 4- V x W - /), t)- < n x ( Vp + Vxw-/),f> 

— (V(V -a;), r)+ < V • n • r > 

+(w — V X u, f) 

—(V x (a; — V x u), v)+ < n x (u> — V x u ), v > 


35 



— (V(V • u), u)+ < V • u, n • v >= 0. (8-2) 

From (8.2) after simplification we obtain the following Euler-Lagrange equations and 
boundary conditions: 

A p = V • / in Cl, (8.3a) 

5 = 0 or n • (Vp + Vxw-/) = 0onr, (8.36) 

Au; — Jjj + V x u = — V x / in ft, (8.3c) 

n x r = 0 or n x (Vp + V xw — f) — OonT, (8. 3d) 

n -t — 0 or V ■ a; — 0 on T, (8.3e) 

An -(-V xw = 0 in Q, (8.3/) 

fixii-O or n x (w - V x u) = 0 on T, (8.3g) 

n-v = 0 or V -u — 0 on T, (8.3 h) 


Equation (8.3) reveals that the least-squares weak formulation corresponds to seven 
second-order elliptic equations and seven boundary conditions in which the original bound- 
ary conditions serve as the essential boundary conditions and some first-order equations 
serve as the natural boundary conditions. In the following we list the combinations of 
boundary conditions for different cases: 

(1) n -u,n XW given 

n ■ ( Vp + V x « - /) = 0, 

V ■ uj = 0, 

n X (u; — V X u) = 0. 

(2) p, n • u, n • <jj given 

n x ( Vp + V x W — f) = 0, 
n x (u; — V x u) = 0. 

(3) p, n x u, n - ZJ given 

n x ( Vp + V x w - /) = 0, 

V • u = 0. 

(4) n x u, n x lo given 

n • ( Vp + V Xu} — /) — 0j 
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V • u; = 0, 

V • u = 0. 


(5) u given 


(6) p, n - a;, n x a; given 


Vp + V xw-/ = 0, 
V • = 0. 


n x (a; — V x u) = 0, 
V * u = 0. 


We emphasize again that the least-squares method based on the first-order velocity- 
pressure-vorticity formulation (5.5) does not need any additional boundary conditions. 
Only if someone would like to use, for example, the finite difference method to solve the 
second-order velocity-pressure- vorticity formulation (8.3a), (8.3c) and (8.3f), should the 
additional natural boundary conditions be included. 

We notice that in the second-order velocity-pressure-vorticity formulation in general 
the solution of the pressure Poisson equation (8.3a) is coupled with the solution of the 
velocity and the vorticity through the boundary conditions. The significant advantages 
of the present second-order formulation are that (1) it guarantees the satisfaction of the 
continuity of velocity and the solenoidality constraint on the vorticity without explicitly 
including these two divergence-free equations; (2) it is suitable not only for the standard 
boundary condition but also for non-standard boundary conditions; (3) the differential 
operator is self-adjoint (symmetrical). 


9. The Div-Curl Method for the Navier-Stokes Equations 

Let us first consider the following first-order system of the Stokes problem with the 
standard boundary condition: 


Vp + V x = / in ffc, 

(9.1a) 

V • u; = 0 in ffc, 

(9.H) 

—uj + V x u = 0 in fl, 

(9.1c) 

V • u = 0 in ft, 

(9.1<f) 

u — uy on I\ 

(9-lc) 
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Of course the boundary data ur should satisfy the global mass conservation: 

J to ■ to r = 0. (9.1 /) 

From Theorem 2 we know that System (9.1) is equivalent to the following system: 


V x ( Vp + V x (jj — /) = 0 in n, 

(9.2a) 

V • (Vp + V x u> — /) = 0 in fl, 

(9.2b) 

to ■ (Vp + V X w - /) = 0 <m r, 

(9.2c) 

V x (w - V x «) = 0 in 9, 

(9. 2d) 

V • (u; — V x to) = 0 in Q, 

(9.2e) 

to x (u; — V x to) = 0 on T, 

(9.2/) 

V ■ ui = 0 in Q, 

(9.2 ff ) 

V • u = 0 in Cl, 

( 9.2k) 

u — ur on T. 

(9.2 i) 

Taking into account V x Vp = 0, V ■ V xu = 0, V-Vxto = 0 and Equality (A.4), System 

(9.2) can be simplified as: 


Au; = —V x / in fl, 

(9.3a) 

A p = V • / in Q, 

(9.36) 

to • (Vp + V x w - /) = 0 on T, 

(9.3c) 

A to + V x ZJ = 0 in SI, 

(9.3d) 

to x (<Z — V x to) = 0 on T, 

(9.3e) 

V • w = 0 in fi, 

(9.3/) 

V • to = 0 in il, 

(9- 3g) 

to = TOr on r, 

(9.3h) 

As explained in §2.3, Equation (9.3p) can be eliminated, since 

Equation (9.3d) and 

(9.3f) and the boundary conditions (9.3e) and (9.3h) guarantee the divergence-free of the 

velocity. 


Since Equation (9.3a) implies that 


A(V • a;) = 0 in 0, 

(9.4a) 
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if we specify that 


V • a? = 0 on r, (9.46) 

then V-w = 0 in fl, that is, the solenoidality of the vorticity vector is guaranteed. Therefore 
we can replace Equation (9.3f) by the boundary condition (9.4b). Furthermore, if we 

replace (9.3e) by that _ 

w — V x w = 0 on T, (9.4e) 

then (9.3e) and (9.4b) all are satisfied. Finally we obtain the conventional second-order 
velocity-pressure- vorticity formulation for the Stokes problems. 


Aa> = —V x / in ft, 

(9.5a) 

A u + V x a; = 0 in fl, 

(9.56) 

w-Vxti = 0onr, 

(9.5c) 

u = ilr on r, 

(9.5d) 

Ap — V ■ / in Cl, 

(9.5e) 

(Vp + V x w - /) = 0 on r. 

(9-5/) 


From (9.5) we understand that the calculation of the velocity and vorticity is decoupled 
from that of the pressure, for this reason in the literature Equation (9.5a) and (9.5b) are 
called the velocity-vorticity formulation. We note that this decoupling does not hold if on 
a part of boundary the pressure is prescribed. 

For the Navier-Stokes equations, Equation (9.1a) in system (9.1) is replaced by the 
following non-linear momentum equation [2 ,15]: 


lo x u + Vp + — V x uj — f in ft, 

Re 


(9.6) 


where p should be understood as the total pressure. Following the same steps as those for 
the Stokes equations we obtain the second-order velocity-vorticity formulation: 

—Aw - V x (w x u) = —V x f in fl, (9.7a) 

Re 


A ?/. V X w — 0 in fl, 
u; — V xu = 0 on T, 
u = ur on r, 
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(9.74) 

(9.7c) 

(9.7<i) 



(9.7e) 

(9.7/) 


Ap + V • (a; x u) = V • / in fi, 


n • (u? x u + Vp + — V x w - /) = 0 on T. 
Re 


10. Conclusions 

The least-squares method based on the first-order differential equations is not only 
a powerful technique for numerical solution, but also a useful tool for theoretical study 
of the div-curl equations and the Navier-Stokes equations. The div-curl equations and 
the Navier-Stokes equations in the first-order velocity-pressure-vorticity formulation are 
not overdetermined. The three-dimensional div-curl equations should have two boundary 
conditions. Since the principal part of the Navier-Stokes equations consists of two div-curl 
systems, four boundary conditions on a fixed boundary are needed (if three conditions are 
given, the dummy <j> = 0 on T should be counted as the fourth one) for three-dimensional 
problems, and two for two-dimensional problems. Four different combinations of non- 
standard boundary conditions are rigorously proved to be permissible for the Navier-Stokes 
problems by using the least-squares method. Consequently, the corresponding least-squares 
finite element method with equal-order interpolations has an optimal rate of convergence 
for all unknowns. The least-squares method and the div-curl method are systematic and 
consistent methods to obtain a high-order derived version of the differential equations 
without generating spurious solutions. Specially, the self-adjoint second-order differential 
equations obtained by the least-squares method automatically satisfy the divergence-free 
equations and are suitable for any boundary conditions. 
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Appendix 


A. Operations on Vectors 


V ■ (qv) = gV • u + Vg ■ r, 

(-4-1) 

V x (gu) = gV x v + Vg x v, 

(-4.2) 

(lx®) = (V x u) ■ v — u • (V xu), 

(4.3) 

V x V x v = V(V • u) — Ar, 

(4.4) 

v x (V x u) — ^ V(u 2 ) — (uV)u. 

Z> 

(4-5) 


B. Green’s Formula 

Assume that u, v and g are smooth enough. Integrating (A.l) and using the Gauss 


divergence theorem lead to 

(V • v, g) + (u, Vg) =< n • v, q > . 

Substituting v — Vp into (B.l) yields 

(Ap, g) + (Vp, Vg) =< n • Vp, g > . (#-2) 

Substituting q = V -u into (B.l) yields 

(V -v, V • u) + (t), V(V ■ u)) =< n - v, V • u > . (-^-3) 

Replacing v by V x v in (B.l) leads to 

(V x v, Vg) =< n • (V x v), g > - (-^-4) 

Integrating (A.3) and using the Gauss divergence theorem lead to 

(V Xu, v) — (u, V x v) =< n x u, v > . (^-5) 

Substituting u = Vg into (B.5) yields 

(V x v , Vg) = - < n x Vg, v >=< Vg, n x v > . (#.6) 
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Replacing v by V x v in (B.5) yields 

(V xu, V x v) — (u, V x V x v) —< n xu, V x v > . (J3.7) 

C. Poincare Inequality 

Let fl be a bounded domain with a piecewise C 1 boundary T, then 


Mi < c{||Vp|U + ( f pdx) 2 } v P e zr’(fi), 

Jn 

(C.l) 

Ml? < C{||V P ||2 + (j pixf) v p€ i'(fl). 

(C.2) 
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